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Starting from the Kadanoff-Baym relativistic transport equation and the multiple 
scattering expansion of the self-energy, we obtain the Boltzmann collision terms for 
any number of participating particles to all orders in perturbation theory within a 
quasi-particle approximation. This work completes a program initiated by Carring- 
ton and Mrowczyhski and developed further by present authors and Weinstock in 
recent literature. 



I. INTRODUCTION 

The study of non-equilibrium phenomena is one of the most interesting problems in 
physics. In particular, the study of the formation and evolution of the quark gluon plasma 
(QGP) in heavy ion collisions is both theoretically and experimentally challenging. The 
current dynamical picture of the evolution of the "fireball" involves many stages. On one 
hand, it is clear from a theoretical point of view that the saturated gluon statefthe so-called 
Color Glass Condensate) describes the initial states of the relativistic nuclei 
On the other hand, it is an experimental fact that the final states only consist of fully 
formed hadrons. In between these two extremes, the hot and dense matter evolves through 
many stages; after the initial collision, the quarks and gluons must scatter enough times 
to thermalize and form a QGP. This lump of QGP then must expand hydrodynamically, 
eventually reaching low enough temperatures to undergo a phase transition to hadronic 
matter. These hadrons interact further until the final kinetic freezout temperature is reached. 
None of these stages are static nor are they in global equilibrium. 

In principle, quantum field theory is capable of explaining all of the above processes: 
Schwinger and Keldysh taught us how to formulate the problem a long time ago P, |^ . Ac- 
tually solving a non-equilibrium problem is, of course, a difficult task. For non-equilibrium 
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problems, the most frequently used theoretical tools are linear response theory, hydrodynam- 
ics and kinetic theory. Linear response theory and hydrodynamics assume that the system 
under consideration is in a near-equilibrium state. In fact these two approaches are related 
in the sense that the tr ansp ort coefficients needed in hydrodynamics can be calculated using 
linear response theory [15|. For these two approaches, their formalism is general and does 
not rely on perturbation theory. Although perturbative calculations are almost always the 
only practical way to calculate the transport coefficients, in principle the validity of these 
approaches depends only on the assumption that the system is near equilibrium. 

Kinetic theory, in this sense, is different. A kinetic theory is a theory of colli ding particles 
and as such, it must rely on the smallness of coupling constants 10|, |l3|, . First of 

all, to have well-defined particle degrees of freedom, collision rates must be small enough so 
that the final state particles of an individual collision have enough time to become on-shell 
(decohere) before being involved in the next collision. Second, the decay rate of the particles 
must be small enough so that the particles become on-shell first and then decay. The strong 
point of kinetic theory is that as long as these conditions are satisfied, one can numerically 
study arbitrarily out-of-equilibrium systems. The success of the so-called 'cascade models' 
or 'transport models' for modelling relativistic heavy ions collisions attests the strength of 
such an approach. 

Current transport models rarely go beyond including 2 — 2 and 2—^3 processes and 
resonance decays. In a recent paper, Carrington and Mrowczyhski addressed the multi- 
particle interactions in kinetic theories starting from the Schwinger-Keldysh formulation of 
quantum field theory (see also Ref . ) . For a QGP and also for the early universe, we 
expect the systems to be dense enough for multi-particle interactions to become important. 
In Ref. for example, it is argued that the thermalization of the fireball in their scenario 
would not be fast enough without gg —>■ ggg scatte ring s. Other studies also stress the 
importance of multi-particle interactions (e.g. 2l[). In the context of transport 

Lscosity i 
0,0- 

The method developped in Ref.|l6[ is therefore important because it provides a system- 
atic way to obtain relativistic transport equations that include multi-particle interactions, 
starting from first principles. Using a g(f)^ + X(j)'^ scalar field theory, these authors demon- 
strated that up to 4-loop self-energies, the collision term in the Kadanoff-Baym equation can 
be re-expressed as the usual gain and loss rates in a Boltzmann equation within a certain 
quasi-particle ansatz. 

As the number of loops increases, the number of relevant diagrams increases very quickly. 
The authors of Ref. [3| managed this problem by developing a symbolic manipulation pro- 
gram to reorganize the diagrams in terms of vacuum scattering processes. The method could 



coefficients, it has been shown that for bulk viscosity in hot scalar theories, number changing 
processes are essential and must be included 
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in principle be applied to any loops and any theories, but at higher orders, the task quickly 
becomes cumbersome. 

The goal of this paper is to generalize the results of Q] using a different method. Instead 

of using a perturbative expansion, we use the multiple scattering exp ansion of the self-energy 

described in Ref . 0] to show that the pattern exhibited in Ref . |ld | is in fact general. 

A summary of our earlier work on this subject has been published in the proceedings 
I I I I 

of Quark Matter 2005 23]. A similar work was also made public j2^ shortly after, whose 

main results included only tree-level scattering matrix elements. The present contribution 

completes these works so that the analysis can be extended to all orders of perturbation 

theory involving any number of participating particles. 

The rest of the paper is organized as follows. To set the notations, and also set the 

stage for the rest of the paper, we present in Section |n] a brief derivation of the Kadanoff- 

Baym equation for scalar field theories, emphasizing the approximations used. In Section 

nil we summarize the derivation of the multiple scattering expansion of the self-energy and 

introduce our quasi-particle ansatz. We present the derivation of the Boltzmann collision 

term in Section IIVI and finally conclude in Section 



II. KADANOFF-BAYM EQUATION 



The goal of this section is to briefly discuss the Kadanoff-Baym transport equation for 
a scalar field theory. Since there are many excellent papers and books on the subject 



16, m 
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271 |28|), we do not attempt another derivation here. Instead we focus 



(e.g. 

on the approximations and assumptions used in the derivation and the interpretations of 
the equation. The notations and discussions in this section largely follow those in Ref. [l^ 
(with a different i convention for propagators, however). The starting point is the following 
Lagrangian density for a real scalar field: 



C(x) 



—m 
2 



2 J,2 



3!^ 



A 

4!' 



(1) 



To discuss the Kadanoff-Baym equation, it is convenient to define the following set of two- 
point Green's functions: 



e{xo - yo) {OH{x)OH{y)) + e{yo - xq) (Oh(i/)0^^(x; 
e{yo - a;o) (6/^(x)6j,(y)) + e{xo - yo) (dH{y)dH{x] 
dH{x)dH{y)) 

dH{y)dH{x))=D>{y,x) 



(2) 
(3) 
(4) 
(5) 
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Here Oh{x) is a Heisenberg picture operator and the average is defined with respect to the 
initial density operator (■ ■ ■ ) = Tr (po ■ ■ ■ )• The superscript c (a) stands for chronological 
(anti-chronological) ordering. In the language of the Schwinger-Keldysh closed time path 
formalism, we have Dq = D^, Dq = D"^, Dp = D"^ and Dq = Dj^, where the super- 
scripts 1 (2) refer to the time branch where the time goes from — oo to +00 (+00 to —00). 
It is also useful to define the retarded and advanced two-point Green's functions: 

tDf{x,y) ^ e{x'-y') ([Oh(x), 6^,(y)]) (6) 
iD^^''{x,y) ^ e{y'-x') ([6^(1/), 6h(x)]\ = tD^S'iv,^) (7) 



In this paper, we will use the symbols G^'^' ■ ■ for propagators and 11^'^' ■ ■ for self-energies. 

The derivation of the Kadanoff-Baym equation starts with the following Schwinger- Dyson 
equation for the contour propagator 



G{x,y) = Go{x,y)+i / d^z / d^z' Go{x,y)U{z, z')G{z' ,y) (8) 

Jc Jc 

where the subscript signifies free field quantities and the time contour runs from —00 to 
00 along the real time axis and then back to —00 slightly below the real time axis. Since 
we want to derive a transport equation, we need to turn this propagator equation into an 
equation for the distribution function /(X, p). To do so, one first needs to turn Eq.® into 
an equation for the Wigner transformed G^'^ . Since the Wigner transform G^'^( X, p ) is 
the quantum analog of the classical distribution function (when it is positive) 0, E^j '^^ 
can make the following ansatz for our distribution function: 

e{p')G<{x,p) = e{p')p{x,p)f{x,p) (9) 

where p{X,p) is the Wigner transformed spectral density 

p{X,p) = y"ci^«e^P" ([0(X + V2),0(X-n/2)]) (10) 

To make further progress, it is necessary to make approximations. More precisely, we require 
weak inhomogeneity of the system with respect to the inverse characteristic momentum (al- 
lowing for the gradient expansion) and the Compton wavelength (justifying the quasi-particle 
approximation) in addition to the weak couplings. These approximations are necessary for 
transport theory, and we will come back to them in more details in Section IIIIDI Using 
these approximations and following the steps in we get the Kadanoff-Baym equation 
for a scalar theory: 

^p(^ + v-v^/(x,p) + v\/(x,p)Vp/(x,p) = e{po)C{x,p) (11) 
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along with the mass-shell condition: 



p'^ -m'^ + V{X,p) = (12) 

where V{X,p) = Ils{X) + Re W^^lXjp) is the Vlasov term (corresponding to an effective 
mass). Here Us is the tadpole self-energy, Ep is the solution of the mass-shell condition (fT^ 
and V = is the velocity of propagation. The collision term is given by: 

c{x,p) = i[n<(x,p)(/(x,p) + i)-n>(x,p)/(x,p)] (13) 

where n^'^(X, p) are the Wightman self-energies. These Wightman functions are "com- 
pletely cut" self-energies and can be interpreted as gain and loss rates. Written in this form, 
the transport equation is similar to the Boltzmann equation of classical kinetic theory. The 
collision term includes a gain term with an appropriate Bose-Einstein enhancement factor 
and a loss term 0]. 

The derivation of the collision term Eq. pH|l in terms of the scattering processes is the 
main goal of this paper. The microscopic gain/loss rates themselves must be computed 
using the vaccum quantum field theory with the in-medium modified mass. The procedure 
followed in Ref. 0| is to directly analyze multi-loop (up to four) self-energy diagrams using 
a symbolic manipulation program and express them in terms of physical processes. In this 
paper, we use another method based on the multiple scattering expansion of the self-energy. 



III. SCATTERING EXPANSION OF THE SELF-ENERGY 
A. Presentation of the Expansion 

The purpose of this section is to present the central ideas involved in the multiple scat- 
tering expansion of the self-energy. We will use it in Section IIVI to express the self-energies 
n<''^(X, p) in terms of scattering amplitudes so that C{X,p) in Eq.(fTT|) can be interpreted 
as a kinetic theory collision term. For general equilibrium theories, the scattering expansion 
was accomplished in Ref.js^]. Here we repeat the central points and explain how it should 
be modified to deal with subtleties associated with medium complications. 

In this section, we use equilibrium quantum field theory to derive the expansion. However, 
as long as we are careful and do not use the KMS condition, the results in this section can be 
readily generalized to the non-equilibrium case by a simple ansatz n{Ep) —>■ f{X, p), where 
/(X, p) is the non-equilibrium phase space density. This is a very simple extension of the 
equilibrium theory to the non-equilibrium case; but as long as the Kadanoff-Baym equation 
is valid, it captures a large amount of physics. 
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The starting point of our derivation is the four propagators^ corresponding to the Green 
functions in the Schwinger-Keldysh contour. Denoting the regular (—00, 00) time branch as 
"1" and the reverse {00, —00) branch as "2" , the free propagators for equihbrium situations 
are given by [2^ : 

G'uip) = G^(p)+r°(p) (14) 
G'22ip) = G^(p)* + r°(p) (15) 
G'M = D'M+r\p) (16) 
GUp) = DU-P) + T'^{P) (17) 

where 

GUp) = , (18) 
Dl{p) = ^(/)27r(5(/ - m^) (19) 

are the vacuum Feynman propagator and the vacuum Cutkosky propagator, respectively. To 
compute the Wightman self-energy using this set of propagators, one must use the following 
Feynman rules 



22 



1. Draw all cut diagrams where the cut separates the two external vertices. Label 11^ (/c) 
all diagrams where the external momentum k enters the unshaded ('1') region. Label 
n^(A;) all diagrams where k enters the shaded ('2') region. 

2. Use the usual Feynman rules for the unshaded side assigning G\i{p) to the uncut lines. 
For the shaded side, use the conjugate Feynman rules assigning 6*22 (p) the uncut 
lines. 

3. If the momentum of a cut line crosses from the unshaded to the shaded region, assign 
Gi2{p)- If the momentum of a cut line crosses from the shaded to the unshaded region, 
assign G^i(p). 

4. Divide by the appropriate symmetry factor. 

To see an explicit example employing the above rules, see Appendix ^ 

Note that each propagator can be separated into a zero temperature part and a finite 
temperature part. Furthermore, the thermal phase space factor 

T\p) = n(|/|)27r5(/-m2) (20) 



^ In this paper, we concentrate on scalar theories because we want to reproduce the results of |l6| . However, 
the multiple scattering expansion is general and does not rely on any specific Lagrangian. 
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is common to all and the zero-temperature part of the propagators are exactly those of 
the vacuum Cutkosky rules [s^. Therefore any diagram in this approach can be naturally 
divided into a thermal part (r°(p)) and a zero-temperature part. 

Using the above two properties of the finite temperature propagators, it is possible to 
expand any diagram as a power series in the number of thermal phase space factor. In this 
expansion, each coefficient of the expansion is a zero temperature scattering diagram [2^ . 

At a first glance, this looks like a fairly complicated expansion. First, for any given 
topology of a diagram, one needs to consider all relevant thermal cuts (i.e. sum of all different 
combinations of 1 and 2 vertices) [2^. Then each of these thermal cut diagrams must be 
expanded again according to the number of F^'s. This expansion becomes manageable only 
if one can show that the coefficients of the same F° factors have a simple meaning in the 
context of the vacuum theory. 

That these coefficients cannot be the usual Feynman diagrams is evident from the appear- 
ance of the cutting rule propagators in Eqs. ()14llf7j) . Hence, the coefficients must correspond 
to the vacuum Cutkosky (cut) diagrams. In vacuum theory, Cutkosky diagrams are needed 
when one needs to work with the unitarity condition 

2\mT = T^T (21) 

where iT = S — 1 is the transition operator. Therefore, roughly speaking 

(i) The part of propagators form Feynman diagrams that correspond to the matrix 
elements of T, whereas the G^* part of G22 propagators form Feynman diagrams that 
correspond to the matrix elements of T'^ . 

(ii) The D^p part of and propagators provide the intermediate states connecting 
T"!" and T matrix elements. 

(iii) The thermal part F'' represents the bath of on-shell particles in the medium that 
participate in the collision process described by the matrix elements of T and 



In Ref. 



22| . it was argued that following this reasoning yields 

X (k, {li}m\T^\{Pj}n){{Pj}n\Tk\k, {li}m) (22) 



The subscripts m and n in {li}m and {pj}n indicate that there are m /j's and npj's in each set. 
In here and the following, we use the label k exclusively for the external momentum entering 
the unshaded side of the self-energy diagram. The subscript k on the operator %. means 
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Included 



Not included 



FIG. 1: An illustration of diagrams that are included and not included in 

{k, {k}m\'T^\{Pj}n){{Pj}n\'Tk\k, {li}„ 

that the external momentum always interacts with the other momentum states. Note that 
we do not write \{{pj}n\%\k, {k}m)\'^ but leave {k, {li}m\'T'^\{Pj}n) and {{pj}n\%\k, {k}m) 
explicitly separated; this is to express that even if the momenta {pj} can be interconnected 
in all possible ways, some diagrams needed to make a squared matrix element are missing. 
The subscript "conn." indicates that the joining of Tf! and Tk does not give disconnected 
pieces, even if by itself 7^ contains both clusters and spectators (see Figs. [T] and [21). A 
cluster is a connected amplitude for a certain subset of initial and final momentum states 
that is completely disconnected from the other amplitudes. Note here that the denomination 
"initial" and "final" when applied to states sandwiching 7^ is slightly different from states 
sandwiching the fully connected scattering amplitude (see Fig. Eland the discussion around 
Eqs. flH^ - fjHHl) ). A spectator is a state that does not interact with other states (see Fig. |21). 
For future reference, we also define "fully connected" as the sum of all diagrams where all 
momentum states are connected together when joining two amplitudes (see Fig. |21 for an 
illustration of this). 



B. Medium Related Complications 

To turn the collision term of the Kadanoff-Baym equation into the usual Boltzmann equa- 
tion collision term, we must determine under what conditions Eq. fl22|) can be approximated 
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FIG. 2: An example of a 3 3 process included in the multiple scattering expansion 
(fcj {^i}2|'^!|{Pj}3)({Pj}3|'^'|^) {^j}2) • The lines cut by the dashed lines are final states and 

conn. 

the arrows indicate the complex conjugated side. The blobs represent connected amplitudes. The 
{pj}3 are the momenta of the final states while the {li}2 are the initial states weighted by thermal 
factors. The external momentum k is always interacting. contains both clusters and spectators. 
The expression T^Tk contains all the possible connections between each other. The (...) represents 
diagrams similar to the ones already shown but with all momenta and the shading reversed. The 
subscript conn, selects all the diagrams above the dashed line. The first diagram of T^Tk is the 
only fully connected diagram. 



by 



ni s\ 

s=0 n=0 



X 



\{{Pj}s\M\k,{k}n)\' (23) 



where |({pj}5|A^|A;, {ii}n)|^ is the square of the fully connected scattering matrix element 
between the initial momenta k and {li}n and the final momenta {pj}s- This is the main goal 
of our paper. Note that M. does not carry the subscript k since it is fully connected. 
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The presence of the medium, however, comphcates this task. Due to the processes oc- 
curing within the medium, the ({Pj}n|'^|{^j}m) piece appearing in Eq. (f2^ does not have to 
be composed of only the fully connected diagrams. Some of the initial momenta {li}m and 
the final momenta {pj}n can interact among themselves without interacting with the main 
piece where k interacts. This is illustrated in Fig. ^ and Fig. |21 Figure ^ also illustrates 
why the factor 



K.T{k,{li}m,{Pj}n) = {k,{li}m\'T'^\{Pj}n){{Pj}n\%\k,{li}m) 



(24) 



cannot be readily interpreted as the square of a fully connected scattering matrix element. 
The right hand side diagram in Fig. ^must appear in the square of the T matrix element, 
yet it is missing in /Ct because these types of diagrams do not appear naturally in the 
self-energy calculation. This fact is taken into account by the "conn." prescription. 

The occurence of 'spectators' and 'clusters' in a T matrix element is well known (for 
instance, see jsil). In vacuum theory, these disconnected pieces do not contribute to the 
self-energy calculation because they represent off-shell vacuum fluctuations which are inde- 
pendent of the scattering process. At finite temperature, however, these in-medium processes 
can no longer be neglected. This is because the disconnected pieces can now represent real 
processes occuring between on-shell particles in the medium. Furthermore, quantum inter- 
ference can occur between these /c-independent in-medium processes and the process where 
k is directly involved. The first diagram in Fig. ^is exactly of this type. 

Another complication is the fact that no momentum state is stable at finite temperature 
since, within a mean free path, the momentum of any particle will likely change. This 

Therefore, 



32 



implies that there are no real asymptotic states at finite temperature 
one must always start with resummed propagators and make an appropriate quasi-particle 
approximation. 

The use of resummed propagators is also essential to deal with so-called "pinch" singu- 
larities. When computing diagrams in finite temperature quantum field theory, one must 
be careful with self-energy insertions since the resulting product of contour propagators 



^for equilibrium sys- 



35l |. However, for out- 



can give rise to ill defined products of delta functions. Fortunate^ 
tems, these singularities cancel due to the KMS condition j22, |32 , 
of-equilibrium systems in which the KMS condition is not satisfied, those singularities do 

at any finite order in perturbation theory. But as pointed 



not cancel 
out in 



34 
37 




, an out-of-equilibrium system is integrated only over a finite volume; 
hus, products of delta functions do not diverge in physical situations. As discussed in Refs. 

38|, the way around this difficulty is to use resummed propagators, where no 



36, 



37 



self-energy insertions (and thus no pinch singularities) appear. To show that Eq. (j23|) is 
indeed a good approximation of Eq. and to deal with pinch singularities correctly, it is 
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thus essential to use resummed propagators. 



C. Resummed Perturbation Theory and the Skeleton Expansion 



The spectral density for a stable particle in vacuum has two separate pieces, an on-shell 
5-function piece and an off-shell continuum piece. Stability requires that these two pieces are 
well separated. In a medium, however, on-shell states are not stable since elastic collisions 
can easily change their momenta. This implies that the zero temperature part and the 
in-medium part may not separate cleanly. 

Unless the external particles are either stable or we assume that the interactions are 
turned off at t = ±00, the scattering matrix elements between on-shell particles are not 
well defined. Technically speaking, the validity of the LSZ formalism depends on the fact 
that we can cleanly amputate the external legs of a Feynman diagram. That is, the self- 
energy insertions along any external leg result only in renormalization of the mass and the 
wavefunction. At finite temperature, this can no longer be cleanly achieved because the 
self-energy has a non-zero imaginary part for any momentum. In other words, both the 
on-shell piece and the continuum piece in the spectral density broaden so that there is no 
longer a gap between them. Nevertheless, if the interaction strength is weak and the medium 
is weakly inhomogeneous, the on-shell piece and the continuum piece should be separated 
enough so that we can use a 'quasi-particle' approximation. We will do so in a short while. 
Here, we first discuss the skeleton expansion whose result is needed later. 

To address the self-energy problem in a systematic way, consider the resummed pertur- 
bation theory where the propagators are fully resummed. In such a theory, the propagators 
in Eqs. f|14llf7|) become: 




Gn{p) 

G22{V) 

Guip) 



Gp{p) + T{p) 
GUp) + r{p) 
0{p')p{p)+np) 



(25) 
(26) 
(27) 
(28) 



where 



Gf{p) 



— rn? + linip) + i\Iii{p)\ 



(29) 



is the resummed analogue of the Feynman propagator and 



pip) 



2n,(p) 



(30) 



Ip2 _ _|_ nij(j9) + in/(j9)p 
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is the resummed spectral density. The thermal factor is given by: 



n{\p |)sgn(p )p{p) 



(31) 



In these expressions, IIr{p) and 11/ (p) are the real and imaginary parts of the retarded 
self-energy and sgn(a;) denotes the sign of x. 

The original multiple scattering expansion of the self-energy (c.f. Eq. (|22|) ) is based on 
an expansion in terms of thermal factors. As mentioned earlier, we need to use resummed 
propagators in order to deal with pinch singularities correctly. Because of this fact, Eq. (j22p 
is not useful in a practical sense. It must be revisited using resummed perturbation theory. 
Moreover, to reproduce "Boltzmann" Wightman's functions (c.f. Eq. ()23|) ). it is convenient 
to modify the original expansion so that thermally weighted initial and final states appear 
more naturally. Thus, we follow a reasoning analogous to the one that led to Eq.()22j) but 
using resummed propagators and organizing the perturbation series differently. In order to 
do that, we rewrite the propagators Gu and G21 as 



The [1 + n{\p^\)]p^{p) factor in G12 and G21 is naturally associated with a final state 
momentum in any given scattering process, and the n{\p^\)p^{p) factor in G12 and G21 is 
associated with an initial state momentum in any given scattering process. The idea is to 
expand the full self-energy using the propagators (jn^ - ljHHI) : see Fig. El for an illustration 
of this procedure. Note that since we use resummed propagators, only skeleton diagrams 
appear in the self-energy. As usual, skeleton diagrams are defined as diagrams without any 
self-energy insertions. The result of this procedure can be written in a very suggestive form: 



Guip) 

G2l{p) 



9{p') [l+n{\p'\)]p^{p)+9{-p')n{\p'\)p^{p) 



(32) 
(33) 



with 




(34) 




X {k, {li}n\Al\{pj}s){{Pj}s\Ak\k, {h}n) 



skeleton 



(35) 



where the resummed phase space factors are given by 



dT^p) 




0{p')p+{pM\p'\) 



(36) 



dr+{p) 



(2vr)4 
d'^p 



^(p>+(p)[l + n(b°|)] 



(37) 



(2vr)4 
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n^(k) = 5^ 



or G21 
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FIG. 3: The first diagram represents the sum over ah skeleton diagrams that have % cut hues. 
The blob on the shaded side is made of G22(p)'s while the blob on the unshaded side is made of 
Gii(p)'s. Cut lines are associated with the cut propagators Gyiip) and G2i(p). The second diagram 
is obtained from the first by expanding the cut propagators using Eq. H38|) and by changing the 
variable so that 9{—p^) — > 0{p^). Thus, all the states weighted by dT~^ becomes final states while 
the states weighted by dT~ becomes initial states. The lines cut by the dashed line are final states. 
The arrows indicate the complex conjugated side. We neglected symmetry factors for simplicity. 



with the relation (true for rotationally invariant systems 



dT+{p) + dT4-p) 



Gl2{p) 



(38) 



Here {{pj}s\-^k\k, {li}n) is purely made of Gu and represents the sum of all amputated skele- 
ton diagrams^ between a set of initial momenta and a set of final momenta {pj}s- The 
factor (A;, {^j}n|^l;|{Pj}s) has the same interpretation and is purely made of G22; however, 
it must be emphasized that 



K,A{k,{li}nAP3}s) = {k, {h}n\Al\{pj}s) {{Pj}s\Ak\k, {li}n) 



(39) 



skeleton 

is not the fully connected matrix element squared that appears in the Boltzmann collision 
term. This will be examined further in the next section. 



D. Approximations 

We know that in order to interpret the right hand side of the Kadanoff-Baym equation 
(c.f. Eq. fll3p ) as the collision integral in the Boltzmann equation, the self-energy 11^ (which 
appears in the loss term) should have the form shown in Eq. ()23|) . Comparing Eq. ()23|) with 
Eq.(j2SI) one can see the following differences: 

^ Note that Ea.lj^ contains no spectator, but still includes clusters. We deal with the latter in Sect. IIII Dl 
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(i) The momentum integrals in Eq. ()23p are over on-shell 3- momenta whereas the momen- 
tum integrals in Eq. (j35p are over 4- momenta. 



(ii) There are some extra terms in Eq. ()39|l . First, JCa is made of the thermal propagators 
Gii(g) and G22{<l) while {{pj}\J^\k, {/j}) contains only the vacuum ones (Cj^i^q) and 
G°*(g)). Also, ICa is made of connected diagrams that contains clusters (but no 
spectators), while {{pj}\J^\k, {li}) contains only the fully connected diagrams (see 
Fig. 121 for the difference between connected and fully connected). 

(iii) There are some missing terms in Eq. ()39p . The terms of the form G*p{q)GF{q) (see 
Fig. and example in Appendix ^ are not part of the skeleton expansion in Eq. (j35|) , 
by definition. These terms were present in the original multiple scattering expansion 
Eq.(j22I), so they appear to be lost in the skeleton expansion of resummed propagators. 
They are actually included in the resummed phase space factor. 

Since the Boltzmann equation is only a classical approximation to the full Kadanoff- 
Baym equation, one should not expect that Hg^j^.^ [k) in Eq.(j2Sl) is exactly the same as 
n^(/i;) in Eq. (jn3|) . Rather, the question to ask is, "In what sense is Eq. ()23|) contained 
in Eq.()35p?" Physically, the Boltzmann equation should be a good approximation in the 
weakly-interacting dilute-gas limit. Hence, it should be no surprise that the answer to the 
above question is provided by a quasi-particle ansatz. 

First consider point (i). The connection between the 4-momentum and the on-shell 3- 
momentum is provided by the following quasi-particle approximation. In a medium, no 
momentum state can be absolutely stable at finite temperature. However, the Boltzmann 
equation is an equation for the on-shell particle density. This means that it is valid only 
when the mean free path is much longer than the range of the interaction, which occurs in 
weakly interacting systems (i.e. the usual approximations of transport theory, see Sect. ITTl 



Kiy m 



and [IQl). This leads to the following condition for the quasi-particle approximation 

- H^j > |Hf I (40) 

where Hj^ means the imaginary part of the thermal self-energy. If Eq . (j4Uj) holds, the on- 
part of the spectral function is well separated from the continuum part (see Fig. |3]and 



shel 



4 li p and it becomes possible to make the following quasi-particle approximation 



p+{k) = 2Iii{\kolk)GF{k)G*p{k) 
^quasi^^^ = 2Ti6{e - m^) + 2Iij{k) Gv{k)G*y{k)9{k'' ^ m%) (41) 
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p(p) 

A 




FIG. 4: Qualitative depiction of the spectral density. In the quasi-particle approximation, the 
left-hand side is approximated by 6{p'^ — m^) while the right-hand side is given by 11/ (p). 

where the 6'-function ensures that we do not double-count the on-shell contribution. Here 
= TV? — Il/j is the thermal mass and we defined 

^^^^^ = p^-ml + rliUp) 

with ny denoting the imaginary part of the vacuum self-energy. Equations (j4H) - (j42j) con- 
stitute our quasi-particle approximation. 

Point (ii) is about clusters in IC^. The diagrams that contain clusters represent interfer- 
ence between the processes that are occuring within the medium and the scattering process 
in which the original external momentum k is involved. For instance, see the diagram on the 
left hand side of Fig. ^ As such, it is not possible to interpret these terms as a matrix ele- 
ment squared because terms representing independent scattering processes (as shown on the 
right hand side of Fig. ^ are missing. It is conceivable that one may introduce these missing 
terms via a multiplicative constant analogous to an effective potential. But in this work, 
we need to ignore these interference terms and note that this is a defect of the Boltzmann 
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equation approach. We simply replace Ka with 



where {{pj}s\M.\k, {li}n) now consists solely of fully connected networks of Gy's 
(c.f. Eq.(I121))- Note that the restriction that no self-energy insertion appear still applies. 
Note also that ignoring these independent scattering processes (equivalent to neglecting the 
thermal factors r(p) in Eqs. (j25|) - (j26|) ) would break gauge invariance in a g au g e theory. As 



is well-known from the Ward identities for Hard Thermal Loops (HTLs) [39|, |40| , both prop- 
agators and vertices must be HTL resummed to preserve gauge invariance. In our case, 
throwing away these independent scattering processes would prevent getting the necessary 
thermal corrections to the vertices needed for preserving gauge invariance. We come back 
to this point and possible solutions in Sect. |Vl 

Point (iii) is concerned with the problem of generating all the missing diagrams in our 
skeleton expansion so that we recover full matrix elements squared. This is a non-trivial 
problem and we devote the next section to its solution. In order to make progress, we need 
to use the following relationships: 

e{p^)n{\p^\)2Iii{p) = e{p^)U<{p) (44) 
9{p^)[l + n(|p°|)]2nj(p) = ^(pO)n>(p) (45) 

Equations pijl - ljl^jl are a rewriting of the usual KMS condition n(p°)n-^(p) = (1 + 
n{p^))Il'^ [p) . This relation is only valid in equilibrium and makes the collision term of 
the Kadanoff-Baym equation vanish identically. However, in an nonequilibrium system, the 
KMS condition is not satisfied and should not be used. We argue here that it is still consis- 
tent to use Eqs. ()44|) - ()45p even if we are in a nonequilibrium setting as long as we only use 
it in Eq. ()46|) for the computation of H^(p) and Il'^{p). The key observation is that H^(p) 
and H^(j9) are proportional to each other up to gradient terms: this is the essence of the 
Kadanoff-Baym equation The effect of these gradient terms would be to add higher 

order gradient terms in our computation of H^(p) and H^(p). These higher order terms can 
be neglected, since the Kadanoff-Baym equation is obtained from an expansion to lowest 
order in the gradients. In this sense, Eqs. (j44j) - (j45j) are approximate expressions that are 
consistent with the the Kadanoff-Baym equation. 

IV. COLLISION TERM DERIVATION 

As explained in Sect. 1111 Dl we need to show how the missing terms in the squared 
scattering matrix elements are generated. With the approximations presented in the last 
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section, Eq.(|H3j) becomes 

n=0 s=0 ' ' k=l ^ ' 



I skeleton 



X (A;,{/fc}„|^^|fe},)({j9,},|^|A;,{4}„)| 

with a similar expression for 11^ (A;). This equation defines a self-consistent equation for 11^ 
and n^. Our problem is now reduced to showing that the Boltzmann ansatz, Eq.()23p. is a 
solution of Eq.(jini)- 

Qualitatively, we can argue as follows: the scattering matrix element {{pj}\M.\k, {h}) is 
now made of zero temperature Feynman propagators Gv{p)- However, we cannot yet say 
that 

/Cm = {kAk}n\M^\{p,}s){{p,}s\M%{l,}^)l^^^^^^^ (47) 

equals {lk}n)^ ■ This is again because diagrams in the skeleton expansion can- 

not include self-energy insertions and hence do not contain all Feynman diagrams con- 
tributing to \{{pj}s\M.\k,{lk}n)\^ ■ We should emphasize here that the matrix element 
{{Pj}s\J^\ki {lk}n) itself does contain all necessary diagrams. It is the joining of the fi- 
nal state momenta {pj}s that makes the distinction between {lk]n)\^ and fCu- 
This is illustrated in Fig. Elfor the part of 11^ which involves 4 particle (2-2 or 1-3) processes. 
The same figure also illustrates that the missing diagrams are the ones with self-energy in- 
sertions. 

Since the right hand side of Eq. pUj) includes 11^ and 11^, this equation must be solved 
iteratively; hence, the missing diagrams must come from this iteration procedure. In other 
words, diagrams missing in the skeleton sum K,m should all be generated by the lower 
order (in the number of legs) skeleton diagrams by trading some on-shell (5-functions with 
Gyll^'^Gy combinations. For example, the diagram labelled (e) in Fig. is generated by 
one iteration as shown in Fig.lHl Fig. [7| contains a more elaborate example. Hence, one can at 
least see that Eq. (j23|) and Eq.(|46|) contain the same set of diagrams. If the symmetry factors 
all work out, then we have shown that the two expressions are equivalent. An example of 
this procedure involving 2 to 2 scattering processes is presented in Appendix 1X1 

Formally, what we need to show is that Eq. lj^Hjl is the iterated solution of Eq. ()46|) . To do 
so, we define 

(ir°(p) = -^0(pO)n(|/|)27r5(p2-m|) (48) 



(27r) 



dVlip) = ^Ap')1^ + n{\p'\)]27r6{p' - m|) (49) 
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FIG. 5: Diagrammatic depiction of I-/W4P. Here C„, denotes the sum of all 1- 
particle irreducible skeleton diagrams with n external legs. Diagrams (a-d) make up 
(^) while (e) is generated by lower order (in the number 

of cut momenta) skeleton diagram as shown in Fig. |SJ 




/ 








/ 
/ 





n 





Iteration 



/ 








/ 
/ 






A part 

that is missing in 
the skeleton sum 



FIG. 6: Lower order (in the number of the cut momenta) skeleton diagram that generate the 
2-particle reducible diagram (e) in Fig. [S] when iterated once. 



and 



dv\v) = dTlip) + rfr° (-P) = ^,gUp) 



and rewrite Eq. (j23p as follows 



00 / „ n 

n=2 ■ V"^ i=l 



(50) 



(51) 
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FIG. 7: A nested non-skeleton diagram that contributes to jMyp. Also shown are the lower order 
skeleton diagram which generates the nested non-skeleton diagram when iterated with M5. 




FIG. 8: Graphical representation of Eo. HSlj) . 

where Mn+i{k, {li}n) is the sum of all Feynman diagrams involving k and a set of n external 
momenta Depending on the sign of the energy, an /j can be either an initial-state 

momentum or a final-state momentum. A graphical representation of Eq. (j3T|) is given in 
Fig. ISl Note here that a similar result has been obtained in the two-loop case 42|. Using 
the symmetry of the bosonic multi-particle state, we also rewrite Eq. ()46|) as 



n=0 s=0 



nl s\ 



00 00 

s=( 

n 

E 

m=0 

s 

E 



nl 



mlin — m) 



t=Q 



t\{s-t)\ 



n (^a) 

a=l 



n 



L(2vr) 



c=l 



n 



\ (2vr)^ 



(52) 



X {k, {la,lb}n\MmPc,Pd}s){{Pc,Pd}s\M\k, 

skeleton 

where the tilde over a momentum variable is a reminder that it is an off-shell momenta. 
Substituting Eq.(|3T| into the right hand side of Eq. (|3^ and shifting some indices, one 
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arrives at 



oo 



oo 



oo 



oo 



n=0 m=0 t=0 s=0 



■ a=l h=l ^ ' 



oo /. "6 




X (/C, {la,lb}n+m\M^\{'Pc,Pd}s+t){{'Pc,'Pd}s+t\M\k, {la,lb}n+m) 



skeleton 



(53) 



What we need to show is that Eq. ()5ip and Eq. ()53p are in fact equivalent. 

To show that the two expressions are indeed equivalent, we start with Eq. (j^ . Consider 
an M„ in Eq. ()5ip . What we would like to do now is to decompose M„ in terms of 1-particle 
irreducible parts. Since Eq. (jHT|l is the result of the skeleton expansion, the lowest order (in 
the number of legs) M„ is M3 and has only the 1-particle irreducible part shown in Fig. |H1 
For n > 3, M„ can be decomposed into the 1-particle irreducible diagram involving all n 
particles (denoted by C„) and diagrams that are combinations of lower order C^'s and Mm''s 
with m < n and m' < n as illustrated in Fig. IHl The Mm' part in turn is made of lower 
order C"s and M's, and the iteration continues until no M is left. 

Since the iteration starts with M3 = C3, one can then say that M„ consists of all tree-like 
diagrams but with the vertices replaced by appropriate Cm with m < n. These diagrams, 
however, are not the tree diagrams of the underlying theory since the tree diagrams in the 
underlying scalar field theory can only have 3-momentum and 4-momentum vertices while 
in these tree-like diagrams any Cm can play the role of a vertex. As an example, we show 
graphical representation of this decomposition for M3, M4 and M5 in Fig. |H1 

The diagrams that contribute to |M„p can be now generated by summing over all pairs 
of the tree-like diagrams from M„ and M* and over all possible distinct joinings of the 
n — 1 legs (excluding k) between the chosen pair; see Fig. El for an example. Among these 
joined diagrams, the majority are skeleton diagrams. A non-skeleton diagram results when 
a tree branch from the M„ side fully connects with a tree branch from the M* side as in 
the diagram (e) of Fig. El These branches could be a primary branch which connects to 
the root C (where the external momentum k enters) by a single line, or an r-th removed 
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FIG. 9: Decomposition of M„ in terms of the 1-particle irreducible part C„ and other Cm and 
Mm with m < n. Here n and m indicate the number of external momenta involved in the process 
including k. For the process involving 3 external momenta, we have M3 = C3. The M4 that 
appears in the third line can be further expanded using the second line. 

branch which needs to go through at least r number of C's to join with the root C . Since 
all possible tree-like diagrams are being summed over, it is clear that these joined branches 
sum to |Mq+ip, where q is the number of joined legs. It is also clear that replacing all joined 
branches with the on-shell part in Eq.()41|) reduces the diagram to a lower order skeleton 
diagram and this reduction is unique. Consequently, 



\MN+i{k,{k}N)\ 
1 




00 00 00 00 

m=0 n=0 t=0 s=0 



e{ll)Gv{h)G*y{h 



6'(Z°) Gv{ld)Gy{ld) 



00 ^ 

J2—^Mu,+l{^b,{hJu,) 



.■"6=2 

00 



.Vd=2 



t\s\ \ J J-l (2tt 
\'' d=i ^ 

X {k, {laJbjn+mlM^ilc, Id} s+t) {{Icjd} s+t\M\k, {lajb}n+m) 

X 6{N — n — s — Ub — Vd) 



M,,+i{idAhJv,)\' 



skeleton 



(54) 



with the external momenta given by {li}N = {la,lc,ku^,lj^Jn+s+ut+va and 6{N - n - s - 
Ub — Vd) is a Kronecker-(5. The symmetry factors in Eq.()54|) are obtained as follows. One 
starts with 

I 

{k, {li}n+m 



(55) 



{n + m)!(s + t)\ 

where the factorial factors accounts for the bosonic symmetry. For the initial momenta 
{li}n+m, there are (n + m)!/(m!n!) number of ways to choose m legs that receives a self- 
energy insertion. Combining with the original l/{n + m)!, this yields the l/(m!n!) factor in 



21 



the second line. The factor in the third hne is similarly obtained. The l/wfe! and 

factors come from the fact that the symmetry factor for a diagram is given by the 



order of the permutation group that leaves the diagram intact 22l \24\. 

Eq.(|54p is exactly the content of Eq. (j53|) . We have thus proved that Eq. (j5ip is the 
iterative solution of Eq.ipUj). 

So far, no equilibrium assumption has entered into our derivation except for Eqs. ()44|) and 
fl45|) . Hence it is possible to replace the bosonic equilibrium distribution function 
by its non-equilibrium counterpart f{X,p). Within our quasi-particle approximation, the 
Boltzmann collision term is given by 

C{X,k) = 

_dH 

n\ s\ 



oo oo 



s=0 n=0 



in + /(A-, k)i E E i ( / n (^/(A-. >.))(/ n 7^11 + /(^^ Pi 



X \ {k,{Pj]s\M\{li}n) 



oo oo 



n\ s\ 

=0 n=0 




X \{{p,}s\M\k,{k}r:)\ (56) 

where |({pj}|A^|/c, is the full scattering matrix squared in vacuum with zero temper- 

ature masses replaced by thermal masses. The sum n + s should, of course, be greater than 
or equal to 2 and the energy-momentum conserving 5-function is absorbed in the definition 
of matrix element squared. This result is similar to the one in Ref. [2^, with the difference 
that in Ref. [2^ only tree diagrams are included in the calculation of |Mp. The analysis in 
Ref. [2J| also did not fully include resummed propagators. 

The main feature of Eq.()56|) is that it is naturally expressed in terms of full scattering 
matrix elements. They can be evaluated theoretically using perturbation theory or by using 
data from experiments. Note also that all the information about the theory is in the matrix 
elements. Therefore, the expressions of the gain and loss rates should be valid for any 
interacting theory of a real scalar field. In the case of cj)^ and 0^, it reproduces the results of 
[lot for 1,2,3,4-loop contributions at tree level. In principle, our analysis could be generalized 
to gauge theories but the preservation of gauge invariance must be carefully dealt with. 



V. DISCUSSIONS AND CONCLUSIONS 

Equations Eq. ()23|) and ()56p are our main results in this paper. Namely, we have shown 
that the Wightman self-energy can be approximately expressed in terms of the thermally 
averaged full scattering matrix elements. We then used this fact to express the Boltzmann 
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are often used in various calculations 
^ and Monte- 



45 



46 



47 



collision term to all orders. Equations Eq. ()23p and 
such as in-medium particle property calculations |4^ 
Carlo simulations of many-body systems; but, as far as we know, a systematic derivation 
has been lacking. Within our quasi-particle approximation, our result is valid for all orders 
of perturbation theory and any number of participating particles. 

Our final results are quite independent of the details of the underlying theory because 
they involve only the scattering matrix elements. Even though our discussions and results 
are for real scalar theories, they can in principle be generalized to more complicated theories 
such as gauge theories. This is clearly the subject of a future work. 

An advantage of having such an all-orders formulation is that these matrix elements do 
not have to be calculated from first principles; instead, they can be measured in vacuum 
and used as an input for the in- medium calculation. Also, having derived Eq. ()56|) . we can 
use it as a basis for Monte-Carlo simulations of many-body systems with confidence. 

While obtaining our main results, we have made some approximations. Here we would 
like to list them and comment on them. 

(i) Approximations in deriving the KadanofT-Baym equation. 

These approximations are valid in the weakly inhomogeneous system and the weak 
interaction strength limit. Since the Kadanoff-Bavm equation has been much discussed 

for some discussion in the context of 



in the literature (see for example 
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2PI effective action methods), we have little more to add to this point. 

(ii) Quasi-particle approximation of the spectral density, Eq. pT|) . 

This is somewhat different than the quasi-particle ansatz made in jioi 24 1 in that we 



have included the continuum part of the spectral density. This approximation should 
be valid under the same conditions that the Kadanoff-Baym equation is valid. To 
go beyond that, we need t o repea t our analysis using the off-shell transport theory 
developed in Refs. jS^ 



ed to re pea t ou: 



53|. This is clearly outside the scope of the 



present paper. 



(iii) Ignoring thermal corrections to the interaction vertices, Eq.( 

Our quasi-particle ansatz could easily accomodate the in-medium mass, but we had 
to make an approximation where we ignored thermal vertex corrections. Physically, 
this amounts to ignoring interference between the purely in-medium processes and 
processes in which the external momentum k is involved. In many physical situations, 
ignoring these corrections can lead to serious inconsistencies. In the case of real scalar 
theories, this is not much of a concern. But when generalizing to gauge theories. 



23 



thermal vertex corrections must be taken into account. For instance, in hard thermal 
loop calculations in hot QCD, it is essential to include thermal corrections to both 
vertices and self-energies to preserve gauge symmetry. 

It should be possible generalize our formulation to include thermal vertex corrections, 
but we will have to pay a price: the scattering matrix elements need to include dis- 
connected pieces so that both diagrams in Fig. [T] are included. Suppose we denote the 
sum of all cut in-medium diagrams (the analogue of the vacuum bubbles) by Wt- In 
vaccum, this is zero because the vacuum cannot spontaneously generate propagating 
particles. In a medium, however, Wt encodes contributions from the independent 
scattering processes occuring in the medium. Just like the vacuum phase, the missing 
diagrams in K-a can be added by multiplying the self-energy by e^^. But this extra 
factor is not a part of the self-energy. Hence, we must compensate it with an additional 
factor of e"'^^. Whether we can repeat our analysis in this way is currently under 
investigation. 



In summary, we have greatly extended our previous work [23[ and a similar work by 



Weinstock 



24| to derive the Boltzmann collision term from the Kadanoff-Baym equation as 



fully as possible. Along the way, we have identified the approximations that have to be made 
and possible remedies to include the ignored effects. Work along this line is continuing. It 
is hoped that this work lays a firm foundation for the inclusion of many-body scattering 
processes in Monte-Carlo simulations. We also mention that it is possible to generalize this 
result to include strong classical sources. First steps in this direction have been taken in 
Ref. 



57| and work along these lines is currently in progress. 
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APPENDIX A: 2 ^ 2 PROCESS EXAMPLE 



In this example, we illustrate how the iteration procedure described in Sect. |TV1 adds 
necessary diagrams to get full scattering matrix elements using the gcfr" scalar theory. We 
will look at the first non-trivial and kinematically allowed process: the 2^2 process. We 
start from the skeleton expansion of the Wightman self-energy. To lowest order, the expan- 
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FIG. 10: Lowest order skeleton diagrams in the g^i^ theory. 

sion includes the diagrams in Fig. ^1 The mathematical expression corresponding to the 
Wightman self-energy is 

"^^^^ = ^/ (0i(0i^2i(j>)G'i2(g)(2vr)V(A: + p-g) 

+ Y / n (|^(2vr)-5^(A: + q,- q,)5\q, - p - q,)5\p + q, - q,) 



X 



G'i2(9i)G'22(?2)G'2i(g3)G'n(g4)G'2i(p) + G'ii(gi)G'i2(g2)G'22(g3)G'2i(g4)G'i2(p) 



+other cuts of 0{g^) + (Al) 

where 'other cuts' do not contribute to the 2 to 2 process. 

Neglecting the thermal phase space factor inside Gu{k) and G22{k) yields 

"^^^^ = f / (^(0i^2ib)G'i2(g)(2vr)^5^(fc + p-g) 

+ y / n ff)4 ^l], {27rY'S\k + q,- q,)5\q, - p - q2)S\p + q-, - q,) 

X Gi2(gi)G'y(g2)G'2i(g3)G'y(g4)G'2i(p) + G'y(gi)G'i2(g2)G'y(g3)G2i(94)G'i2(p) 

(A2) 



where we defined 



— rriT + I 



Let's consider the first term. 

,2 r ^4^ 



= 1/ (2^4 ^^J^, G,,{p)G,M{2^mk + p- q) (A4) 
Using the fact that the spectral density is odd {p{—p) = —p{p)) and by making a suitable 
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change of variable, this becomes 

X { [1 + n(|/|)] [1 + n(|g°|)] 5\k - p - q) 
+ [l+n(|/|)] n(|g°|)5^(A:-p + g) 
+ n(b°|) [l+ni\q'^\)]5\k + p-q) 
+ + (A5) 

Each terms corresponds to a lowest order transition amplitude between off-shell states com- 
ing from the skeleton expansion as defined previously in /C^. 
The second term is 

nr2)W = n (fj)4 (2l^4 i^^yS'ik + q,- q,)S\q, -p- q,)5\p + q,- q,) 

X [GMG*v(Q2)GMGy(q,)GM + Gv(qi)GMG*y(q3)GMGM] 

(A6) 

We can use the properties of the spectral function and the explicit expression of the cut 
propagators to write the preceding equation as 

W2)W = -W / W ^ {k + qi-q2-q3) 



(2)1-; - Y J (2^(2^(2^ 

xe{q',)e{q'M)p{Qi)p{Q2)p{qs)nM) + [l + ^dg" 

X G*y{q3 - qi)Gv{q2 - qi) + G^(?3 + q2)Gv{q2 - qi) + Gyiqs - qi)Gv{q2 + qs) 

+ G*y{q2 - qi)Gv{q3 - qi) + G*y{q2 - qi)Gv{q2 + q^) + ^^(^3 + q2)Gviq3 - qi 
+ (other processes) (A7) 

Each term corresponds to a different channel of a 2 2 scattering process. The neglected 
terms lead to other processes like 1 <-> 3. If we now make the lowest order quasi-particle 
ansatz, p{p) — > (27r)5(p^ — mf,), the above expression can be written as 

1 f d^qi d^q2 d^qs ,^ .4.4, 



xn{E,) [l+n{E2)] [1+^(^3)] 

X M;Mu + M*Mu + M;Ms + M*Mt + M*Ms + M*Mt 

-|- (other processes) (A8) 

Here M^, Mj, represent the scattering amplitude in the s, t.u channel, respectively. The 
terms in the square bracket is not yet the full scattering element squared since the square 
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terms are missing. We only have the interference terms. The squared terms should come 
from iterating Eq. ()A5|) . We can expand the spectral density for this one-loop diagram using 
Eq. ipTj) and keep only the relevant terms for 2 ^ 2 process at lowest order. This yields, in 
addition to the pure one-loop contribution, 

X [2{2'k)5{p^ - ml)G*y{q)IiM)Gv{q) + 2{2Tx)5{q^ - ml)G*y{p)Yii{p)Gv{p)\ 

X I [1 + n{\p% [I + n(|g°|)] 5\k - p - g) + 2 [l + n(|/|)] n(|g°|)5^(A; -p + q) 

+ n{\p^\)n{\q''\)5\k + p + q)^ 
-|- (other processes) (A9) 

Upon using Eqs. ljl^ and (jl^ this becomes 

X { [1 + n(|/|)] (2vr)5(p2 - ml)G*y{q)Ii> {q)Gy{q)5\k - p - q) 
+ [l+n(|/|)] {2'n)6{p'-ml)GUq)^''{q)Gv{q)6\k-p + q) 
+ n(|p°|)(27r)5(/ - ml)G*y{q)Ii> {q)Gv{q)5\k - q + p) 
+ n(|pO|)(27r)5(/ _ ml)G*yq)^<{q)Gv{q)5\k + p + q)] 

+ (other processes) (AlO) 

The last equation is a self-consistent equation for the self-energy. To solve it, we can use an 
iteration procedure. Since we are considering only the 2^2 processes to the lowest order, 
this means in the right hand side we should replace 11^'^ (g) with the purely one-loop result 

nf,_,)(fc) = (21^4 ^^^y e{p')e{q'){2irf5{p^ - ml)5{q^ - m|) 

x{ [l+n{\p'\)\ [l+n{\q% 5\k - p - q) 
+ 2[l + n(|/|)] n{\q''\)5\k~p + q) 

+ n(|p°|)n(|g°|)5^(A; + p + g)} (All) 

There is an analogous equation for n^_-^^(A;). Upon substitution, Eq. ()A10|) can be shown to 
contain the missing square terms 

1 f (fqi d^q2 (fqs x4.4. 



= 2 J (2.)32ij;,(27r)32E,(2vr)32i^3 ^'"^ ^ (^ + ^1-^2-^3) 
X n{E^) [1 + n{E2)] [1 + ^1(^3)] [m^M^ + M^M, + M^Mt 
+ (other processes) (A12) 
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If we combine our two answers, we get 



u>{k) = n>)(A;) + n>)(fc) 

1 f d?qi d^q2 



2 J (27r)32Ei (27r)32E2 (271)32^. 

Xn(Ei) [1 + ^(^2)] [1 + ^(^3)] 



{2'KY6\k + q,~q2-q,) 



Mu + Ms + Mt 



(A13) 



This is in accordance with our general result expressed in equation 
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